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Abstract 


We apply the rational radial basis functions (RRBFs) method to solve 
the Allen—Cahn (A.C) equation, particularly when the equation has a so- 
lution with steep front or sharp gradients. We approximate the spatial 
derivatives by the RRBFs method. Then we apply an explicit, fourth- 
order Runge-Kutta method to advance the resulting semi-discrete system 
in time. It is well known that the A.C equation has a nonlinear stability 
feature, meaning that the free-energy functional is reduced by time. The 
presented method maintains the total energy reduction property of the A.C 
equation. In the end, five examples to confirm the efficiency and accuracy 
of the proposed method are provided. 
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1 Introduction 


In this paper, we examine the Allen—Cahn (A.C) equation as follows: 

WwW, = Vw, — h(w), z € [a,b], t € (0,T], (1) 
with the initial condition 

w(z,0)=g(z), 2z€ [a,b], (2) 
and Dirichlet boundary conditions as 

w(a,t)=l(t), w(b,t)=k(t), t € [0,T], (3) 


where the parameter v is a positive constant and usually represents the in- 
terfacial width. The term h(w) = H (w) with H(w) = 1(w? — 1)? is a given 
energy potential. 

The A.C equation was initially given by Allen and Cahn [1] to illustrate 
the motion of anti-phase boundaries in crystal solids. The A.C equation is 
a specific type of nonlinear partial differential equation used extensively to 
model various phenomena, such as crystal growth [3, 12], phase transitions 
[6, 28], image analysis [2, 8, 20, 26], texture growth [9, 23, 24, 27], and surface 
dynamics in material science [13, 21]. 

We know that the solution w(z,t) to the A.C equation has the feature 
that the total energy E(w) is reduced over time, and it is given as follows: 


Er(w) = f? (Avlw.? + H(w)) dz, (4) 


and may be used to assess the accuracy of the rational radial basis functions 
(RRBFs) method presented here. 

In general, the analytical solution of the A.C equation is not available. So 
many researchers tend to use numerical methods, such as finite difference [22], 
finite element [11], Fourier spectral with periodic boundary condition [25], 
and radial basis functions (RBFs) methods [29]. In recent years, researchers 
have been interested in developing numerical methods that can preserve the 
energy dissipation property of the A.C equations. 

Among these methods, the meshless technique avoids mesh generation. It 
uses scattered points instead of meshing the domain of the problem. For more 
description, see [14, 15, 16, 30, 36, 38]. The RBFs methods are substantial 
instruments for scattered points interpolation and solving partial differential 
equations (PDEs). When the underlying function or the solution of PDE is 
sufficiently smooth, RBFs methods can produce exponential accuracy, but if 
the underlying function or the solution of PDE has steep gradients or discon- 
tinuities, then the RBFs method may produce solutions with oscillation. In 
such situations, the RRBF'’s method can be used to approximate derivatives of 
functions and to solve PDEs [35]. RRBFs method approximating solutions in 
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nonlinear spaces generated by RBFs is more computationally expensive than 
the standard RBFs methods. For smooth problems, the accuracy obtained 
from these methods may not be worth the additional computations. How- 
ever, for problems with discontinuities and steep fronts, these methods will be 
considerably more accurate, and the additional computations will be justified. 


In 2009, the RRBFs method was used to interpolate functions with steep 
gradients by Jakobsson, Andersson, and Edelvik [17]. Sarra and Bai [35] ex- 
tended the method of Jakobsson et al. to interpolate discontinuous functions 
and solve Burger’s equation. The RRBFs-based partition of unity interpola- 
tion was used by De Marchi, Martinez, and Perracchione [4], and in [37], the 
authors used the RRBF's method for solving the Sine-Gordon equation. 

In this paper, we apply the RRBFs method for solving (1) with conditions 
given in (2)-(3), particularly when the equation has a solution with steep 
front or sharp gradients. We approximate the spatial derivatives by the 
RRBFs method. Then we apply an explicit, fourth-order Runge-Kutta (RK4) 
method to advance the resulting semi-discrete system in time. 


This paper is organized as follows: In Section 2, we introduce the RBFs 
and RRBFs methods. The method implementation is given in Section 3. We 
analyze stability issues in Section 4. In Section 5, we present the results of 
numerical examples. Section 6 is dedicated to the conclusion. 


2 RBFs and Rational RBFs interpolation 


2.1 RBFs interpolation 


Let Tr C R™ be a bounded set, let Z = {z{,..., 2} CT be aset of N distinct 
points, hereinafter referred to as centers, and let F = {f(z{),..., f(zy)} be 
a set of function values. An RBF 


W(z) =v (lz — 2% lla, ae z,2°€R”, 


is a function of one variable r = ||z— z°||2 that is centered at z°, which € is a 
free parameter and called the shape parameter [10, 35]. Also, in this paper, 
m= tl. 


The inverse quadratic RBF 


1 


V(r) = 1 (er)? 


is a strictly positive definite RBF that we use it in all examples. 
The RBF interpolant takes the form 
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N 
p(z) = y cv (llz— Z5ll2,€), (5) 
j= 
which the coefficients c; are obtained by solving the linear system 
Ac =f, 
based on the interpolation conditions p (z$) = f;, where f = [f(z°),..., f(2))]’- 


The elements of the matrix A are of the form 
aig = (l2F- 22,6), 49=1,...,N. (6) 
Moreover, A is a symmetric positive definite matrix and thus invertible. 
The evaluation of the interpolant (5) at M points z; is done by multiplying 


c by E, where the elements of the matrix E are of the form 


eg = (l2i—Zllo,€), t=1,0.,M, 9 =1,2,...,N. (7) 


The first and second derivatives of RBF interpolant are of the form 


-Yop v (llz—- lle, € €)). 
Thus 
N 
2§)) = > /¢; Dy (lef — 2lle, €), 
j=l 


that is, 
Df~ Epe, 


where the elements of Ep are of the form 


(Ep),; = Dy (lei — 2lle, €); i,j =1,...,N, 


and 
-Yop (Dy (ll — 2§ll2, €))- 
Thus 
D(D -Yop (Dy (ll2f - lla, €)), 
that is, 


D(Df) = Eppe, 


where the elements of Epp are of the form 
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(Epp); = D (Dv (lz — 52,6), 49 =1,...,N. 


The following theorem [4] gives an error bound in terms of the power 
function Py,z, where Ni,(I) is the native reproducing kernel Hilbert space 
corresponding to the symmetric positive definite RBF [10]. 


Theorem 1. Let w € C([ xT) be a strictly positive definite RBF, and 
suppose that Z = {zf,...,z29} CT is a set of distinct points and that p is 
the interpolant to f € Ny(L). Then for all z € I, we have 


If (2) — p(2)| S$ Py. 22ers: 


2.2 Rational RBFs interpolation 


The RBFs interpolants have problems, such as ill-conditioning, especially 
when the shape parameter tends to zero. These problems might lead to 
inaccurate solutions, particularly when the functions with steep front or sharp 
gradients are considered [33]. Moreover, when the functions have steep fronts 
or sharp gradients, the rational RBFs method approximates them better than 
the standard RBFs method [35]. These are the main reasons that we use 
rational RBFs interpolants. 
The RRBF interpolant of function f is given by 


and u(z) and v(z) are the RBF interpolants 


N 
u(z) = Soc! o(llz- Zl, 6), 
j=l 


and 
N 


vo(z) = Soc}  (lz—-25lle, €)- 


j=l 


Let u = [u(z),...,u(z%)]’ and v = [v(z?),...,v(2%)|” . By applying the 
interpolation conditions, we have a system of equations that is underdeter- 
mined. Thus in order for the rational interpolant to be defined uniquely, we 
add an additional condition (for more descriptions, see [17, 35]), which leads 
the native space semi-norms [10] of the RBF interpolants u(z) and v(z) to be 
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minimized. By adding the condition, we will have a minimization problem 
with the solution v, which is the eigenvector corresponding to the smallest 
eigenvalue problem 

Rv = »v, 


R= diag (y (Gat) (Fata), (8) 


and A is the RBF system matrix given in (6), f = [f(z{),.. fay, and 
D = diag (f(z{),.--, f(2%)). Moreover, division is elementwise and f? is an 
elementwise squaring of the elements of the vector f. When v is found, u is 
obtained by u = Dv. When u and v are found, we obtain the coefficients of 
the RBF interpolants by solving two linear systems, 


where 


Ac“ =u and Ac” =v. (9) 
Now we evaluate the rational interpolant at M points z; by 


Ect 


= 10 
q- =, (10) 
where division is elementwise, E is the RBF evaluation matrix (7), and Q = 
[Q(z1),..-, Q(zm)]”. 


Now, we calculate the derivatives of the rational interpolant at N points 
z° by applying quotient rule as below: 


(Ac’) . (Epe") — (Ac") . (Epe”) 


Pe 11 
Q a (11) 
Q’ _ 2(Ac").(Epe”)? + (Ac’)?.(Eppc") 

7 (Ac?) 
(Ac").(2 (Epc") .(Epe”) + (Act). (Epve’) ) 
(12) 


(Ac’)? 
F , , tT 37 ; : [ 
where Q = [2 (zf),...,2 (2%)] ,Q= [2 (zf),...,2 (2%)] , Ais the 
RBF system matrix, and Ep and Epp are the first and second derivatives 
of evaluation matrix at N points z¢. 

It is concluded from (9) and (10) that the RRBFs interpolant is made by 
the RBFs interpolation matrix A. Therefore, we are able to provide error 
bounds for the RRBFs interpolation [4]. For this purpose, we must consider 
u and v as the values obtained by sampling some functions r and s € N,(T), 
respectively. 


Proposition 1. Suppose that T C R™, w € C(I xT) is a strictly positive 
definite RBF, and let Z = {z{, z§,...,zQ} CT be a set of distinct points 
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and let Q be the RRBF interpolant to f € Ny(I). Moreover, let us suppose 
r,s €N,(I). Then for all z € I’, we have 


|f(z) — Q(z)| < wa (F(Z) I slau + Irv) Paz (2)- 


So Proposition 1 gives the error bound for the RRBFs interpolant in terms 
of data values and power function, as for the standard RBFs interpolant. 


3 RRBFs method for solving A.C equation 


Now, we use the RRBFs method to solve the A.C equation. In fact, we de- 
velop RRBFs to study the A.C equation with Dirichlet boundary conditions. 
In our method, we employ the derivatives of the RRBFs to approximate the 
spatial derivatives. Then we apply an explicit, RK4 method to advance the 
resulting semi-discrete system in time. In fact, we first approximate the 
spatial derivatives at N centers z¢ at time t”, that is, 


n Cc n Cc n Cc n T 
Wrz (2,¢ )= [Wz (z{,t ), We (25,t \yeedg Woe (zy, t )] ’ 


of the A.C equation 
Wt = VWzz — h(w) 
by the RRBFs method, which is given by computing (12) when 


n c yn c yn ec yn\ lt 
WwW = [W (z{,t ), W (25, ¢ ),..., W (zy, t )] ry 


is considered instead of f. This leads to a system of ordinary differential 
equations of the form 
W; = G(W,t). 


Then we apply RK4 method to solve the above system of equations as follows: 
K, = AtG(W”,t”), 


Ky = AtG(W” + 0.5k1, t” + 0.5At), 


K3 = AtG(W” + 0.5ko, t” + 0.5At), (13) 


K, = AtG(W" + kg, t” + At), 


wetta~wre g(Ki + 2K2 + 2K3 + Ky), 


where t” = t”~! + At and At is the time step size. Clearly, to obtain W! 
initial value, W° is required in the iterative computation. Using the initial 
condition in (2), we have W° = [g(z{), 9(z§),--. (2Z%)] 
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To apply the Dirichlet boundary conditions (3), we directly replace 
W (z{,t”) with w(a,t”) = I(t”) and W (2, t”) with w(b, t”) = k(t”). 


A Stability issue 


Since the RRBFs interpolant is constructed by using the standard RBFs inter- 
polation matrix, the RRBFs method can also suffer from instability problems, 
especially for « + 0. In fact, serious problems of ill-conditioning may occur by 
choosing the wrong values of the shape parameter, particularly for infinitely 
smooth RBFs. In order for the system matrix B to be well-conditioned, the 
shape parameter € must not be too small, but small shape parameters are 
required to obtain better accuracy for the interpolation with RBFs. This 
is known as the uncertainty principle. We choose the shape parameter € so 
that the collocation matrix B has a condition number, «(B), in the range 
10'8 < «(B) < 10", in order to determine the proper value for ¢ [34]. These 
bounds for the condition number are valid when using a computer that im- 
plements double precision floating point arithmetic, but the bounds will be 
different when using other floating point number systems. 


5 Numerical results 


In this section, we consider five numerical examples of the A.C equation to 
validate the presented scheme. The accuracy and efficiency of the method 
are confirmed by calculating D.., L2, and RMS error norms as follows. We 
also investigate energy dissipation property in our numerical experiments to 
illustrate the efficiency of the presented method. 


Loo = ||W— wl.., 
Ly = ||W— wll, 


RMS = / wllW — wll, 


where w and W are the exact and numerical solutions, respectively. Moreover, 
in all numerical examples, we have used uniformly spaced centers 


Also in the examples, we evaluate the total energy Fr given in (4) using the 
composite trapezoidal rule for integration. 


Example 1. In this example, we examine the A.C equation with three large 
gradients 
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Wt = VWzz tw — Ww, z€[-1,1], 


and with the initial and Dirichlet boundary conditions [18] 
w(z,0) = 0.52 —0.5sin(1.57z), w(1,t) = —w(-1,t) =1. 


The results are obtained using vy = 10-°, k = At = 0.01, shape parameter 
€ = 6, and N = 100 centers. In Figure 1, we present the approximate solution 
at time t = 10 and the space-time graph of the approximate solution. Table 
1 shows the values of the energy at various times. It can be seen from Table 
1 that energy reduces over time and that the RRBFs method resolves the 
problem accurately. 


Table 1: Energy of the A.C equation for Example 1 at several times. 


Time (t) t=0 t=1 t=2 t=5 ¢t=10 t = 20 
Energy (Er) | 0.3403 0.1941 0.0738 0.0017 7.6908e-04 _7.6900e-04 


Example 2. In this example, we examine the A.C equation 


Wt = VWzz + WwW — w, 
in the computational domain [—1, 1] x [0,10] with the parameter v = 10°. 
The initial condition is given by [31, 5] 


w(z,0) = 0.62 + 0.4sin (Ge 5p ) ; 
and the Dirichlet boundary conditions are 
w(1,t) = —w(-1,#) =1. 


We solve the problem with the RRBFs method with k = At = 0.01, N = 100 
centers, and the shape parameter ¢€ = 6. We have obtained the solution up 
to time t = 10. The graph of approximate solutions at times t = 0,2, 4,6, 10 
and space-time graph for ¢ € [4,10] are shown in Figure 2. As it is seen, the 
solution has a steep front at different times. Table 2 shows the energy at 
several times. It can be seen from Table 2 that energy reduces over time as 
theory guarantees and that the RRBFs method resolves the problem accu- 
rately. 


Example 3. Consider the A.C equation 


We =VW2zz2 +w—w, z€[0,1], t>0 
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A 0.8 -0.6 -04 -0.2 0 0.2 04 0.6 0.8 1 
Zz 


(a) RRBFs solution at time t = 10. 


(b) Space-time graph. 


Figure 1: Numerical results for the solution of Example 1. 
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(a) RRBFs solution at time t=0. (b) RRBFs solution at time t = 2. 


2 2 
15 15 
1 1 
0.5 0.5 
Xo Soo 
= = 
05 05 
1 1 
1.5 1.5 
2 2 
1 08 06 -04 02 0 O02 04 06 O08 1 1 08 06 04 02 0 O02 04 06 O08 1 
z: 2: 


(c) RRBFs solution at time t=4. (d) RRBFs solution at time t = 6. 


? 
Hi 


. P | 
i Ny i 

; ait Hl 

15 wu 0.5 , 

6 ~ 0 
-0.5 

"4 08 06 04 -02 oO 02 0.4 06 08 1 t 7 1 z 

(e) RRBFs solution at time t = 10. (f) Space-time graph. 


Figure 2: Numerical results for the solution of Example 2. 
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Table 2: Energy of the A.C equation for Example 2 at several times. 


Time (t) t=0 t=2 t=4 t=6 t=10 
Energy(Er) | 0.3161 0.0945 0.0154 5.3755e-04 _7.1318c-05 


with the initial and Dirichlet boundary conditions [7] 
w(z,0) =sin(77z), w(0,t) = w(1,t) =0. 


For the numerical solution, we use N = 100 centers, v = 0.0025, k = At = 
0.001, and the shape parameter ¢€ = 8. The plots of the numerical solution 
and the energy are shown in Figure 3. It can be seen that energy reduces 
over time, as the theory guarantees. Compared with [7], better accuracy is 
obtained. 


Example 4. Consider the A.C equation 


Wt = Wzz + —a (ww), 


in the region 0 < z < 15, where the parameter v = 0.08, and the initial 
condition is given by 


z—-2 
,0) = 0.5 — 0.5 tanh : 
oe) ‘ (=) 


The analytical solution is given in [19] as 


w (z, t) = 0.5 — 0.5 tanh (=) ‘ 
2/2v 

where r = we we obtain the Dirichlet boundary conditions from the exact 
solution and solve this problem with the RRBFs method with k = At = 
0.00001, NV = 200 centers and the shape parameter ¢€ = 2.25. 

The graphs of exact and approximate solutions at time t = 0.04 are shown 
in Figure 4. As it is seen, the solution has a steep front. Table 3 shows the 
Lig, L2, RMS errors, and energy 


Pr(w) = [ (Slw-P pale’ 1)*) dz 


at several times. It can be seen from Table 3 that the energy reduces over 
time and that the RRBFs method resolves the problem accurately. Also, the 
results of the method are well matched with analytical solutions. 
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(c) RRBFs solution at time t = 10. 


(e) Space-time graph. 
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wet) 


wet) 
ee oe 


(d) RRBEFs solution at time t = 40. 


(f) energy evolution. 


Figure 3: Numerical results for the solution of Example 3. 
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Table 3: Errors and energy of Example 4. 


t Lig Ly RMS Energy (Er) 
0 0 0 0 509.2856 
0.001 | 2.3112e-05 2.3259e-05 1.6447e-06 508.2498 
0.005 | 5.8036e-04 1.2000e-03 8.2960e-05 504.0963 
0.01 | 5.8573e-04 1.2000e-03 8.2963e-05 498.9173 
0.04 | 7.8661e-05 1.6396e-04 1.1594e-05 467.8536 
3 rational i 
09 | act 7 
08 | : 
07 | 
| 5 
a0 ‘ : 
gos | ‘7 
04 , 
| 3 
03 + ; 
0.2 1 
1 cpeenemtaidiemataiem:. de ] 


(a) RRBFs and analytical solutions at 


t = 0.04. 


(b) Pointwise error at t = 0.04. 


Figure 4: Numerical results for the solution of Example 4. 


Example 5. Consider the A.C equation 
Wt = Wzz + > 
V 


with a forcing term h added by imposing analytical solution w(z,t) = 
e'cos(27z) in the domain 0 < z < 1, where the parameter v = 0.1 [32]. 
We obtain the Dirichlet boundary conditions from the analytical solution 
and solve this problem with the RRBFs method with k = At = 10~*/8, 
N = 100 centers and the shape parameter ¢ = 8. The graphs of exact and 
approximate solutions at time t = 0.1 are shown in Figure 5. Table 4 shows 
the Do, L2, RMS errors, and energy at several times. It can be seen from Ta- 
ble 4 that the energy reduces over time and that the RRBF'’s method resolves 
the problem accurately. Also, the results of the method are well matched 


2 


(w _ w*), 


with analytical solutions and somehow better than the earlier works. 
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Table 4: Errors and energy of Example 5. 


RRBF RBF [32] 
1 Le fs RMS fies Les Ls 
0 0 0 0 19.2446 0 : 


0.001 | 1.6671e-07 2.6333e-07 2.6333e-08 19.2374 | 2.1880e-06 - 

0.005 | 1.7437e-07 3.7249e-07 3.7249e-08 19.2095 | 2.3599e-06 - 
0.01 | 1.7140e-07 4.2040e-07 4.2040c-08 19.1766 | 2.3762e-06 | 8.3500e-06 
0.1 1.1485e-05 8.0532e-05 8.0532e-06 18.8965 | 1.2587e-05 | 6.6300e-05 
0.3 | 4.2165e-05 2.5234e-04 2.5234e-05 19.5199 | 6.551 1e-04 - 


wat) 


0 01 02 03 04 O05 06 O7 O08 og 1 0 01 02 03 04 O05 06 O7 O08 o9 1 
z Zz 


(a) RRBFs and analytical solutions at —_(b) Pointwise error at t = 0.1. 
t=0.1. 


Figure 5: Numerical results for the solution of Example 5. 


6 Conclusion 


This paper investigated the application of the RRBFs method to find the nu- 
merical solution of the A.C equation with Dirichlets’s boundary conditions, 
especially when the equation has a solution with steep front or sharp gradi- 
ents. The combination of RRBFs in space and RK4 in time gave accurate 
and reliable solutions to the A.C equation. To illustrate the accuracy and ef- 
ficiency of the presented method, we carried out five examples and presented 
the corresponding graphs and tables, and compared the results with some 
earlier works. Moreover, it can be observed that in the case where the an- 
alytical solution exists, the numerical solutions described the good accuracy 
of the method. These results support the assurance of using the method for 
equation (1) when the analytical solution is not known. In addition, for all 
examples, the energy dissipation law was investigated. In order to generalize 
the RRBFs method implementation to higher dimensional time-dependent 
PDEs, we need to localize the method with a modified partition of the unity 
method. It is behind the scope of our paper, and we leave it to our further 
work. 
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